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We present a study of the topological susceptibility in lattice QCD with two degenerate flavors 
of dynamical quarks. The topological charge is measured on gauge configurations generated with a 
renormalization group improved gauge action and a mean field improved clover quark action at three 
values of /3 = Q/g^ , corresponding to lattice spacings of a ~ 0.22, 0.16 and 0.11 fm, with four sea 
quark masses at each (3. The study is supplemented by simulations of pure SU(3) gauge theory with 
the same gauge action at 5 values of (3 with lattice spacings 0.09 fm < a < 0.27 fm. We employ a field 
theoretic definition of the topological charge together with cooling. For the topological susceptibility 
in the continuum limit of pure SU(3) gauge theory we obtain Xt^* = 197ljg MeV where the error 
shows statistical and systematic ones added in quadrature. In full QCD xt at heavy sea quark 
masses is consistent with that of pure SU(3) gauge theory. A decrease of x* toward light quark 
masses, as predicted by the anomalous Ward-Takahashi identity for U(l) chiral symmetry, becomes 
clearer for smaller lattice spacings. The cross-over in the behavior of Xt from heavy to light sea 
quark masses is discussed. 
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I. INTRODUCTION 



The topological structure of gauge field fluctuations, in particular instantons, has been invoked to explain several 
important low energy properties of QCD including the breaking of axial U(l) symmetry and the large mass of the 
J]' meson. Numerical simulations on a space-time lattice provide a non-perturbative tool for the study of these 
phenomena beyond semiclassical approximations. 

Lattice studies of the topological susceptibility Xt as a measure of these fluctuations have been mostly carried out 
for pure gauge theory without the presence of dynamical fermions |0| . Recent determinations by various groups using 

different methods have led to a consistent value in SU(3) gauge theory of xj^ — 200 ± 18 MeV 

Sea quark effects on the topological susceptibility have been much less studied, although dynamical quarks are 
expected to have a strong influence on xt leading to a complete suppression for massless quarks. From the anomalous 
Ward-Takahashi identity for U(l) chiral symmetry, the topological susceptibility is predicted ||-||] to obey for small 
quark masses in the chirally broken phase, 

5] 

Xt = — + O(m^), E = - lim lim (V'V')- (1) 

Pi j m,^OV^oo 

It is an interesting question to investigate whether lattice data confirm a suppression consistent with Eq. (0). 

Pioneering attempts to calculate xt in full QCD were restricted to small statistics and were plagued by long 
autocorrelation times. Progress in the simulation of full QCD, as well as increase of available computer power in 
recent years, has enabled this question to be readdressed with a higher accuracy. A number of pieces of work have 
been reported recently J@-^ coming to different conclusions whether the topological susceptibility is consistent with 
the prediction of Eq. (Q). A common shortcoming in Refs. [p|~p^ is that they have been made at only one lattice 
spacing. Ref. ||l3|| , on the other hand, used only one bare quark mass aruq at each coupling constant /3. 

In this article we attempt to improve on this status by calculating the topological susceptibility in full QCD with two 
flavors of dynamical quarks at four sea quark masses at each of three gauge couplings. We perform calculations on con- 
figurations of the CP-PACS full QCD project |l^. These have been generated on the CP-PACS parallel computer |l^] 
using a renormalization group (RG) improved gauge action |]l^ and a mean field improved Sheikholeslami-Wohlert 
clover quark action |lj]. The efficacy of this choice of action over the standard action has been demonstrated in 
Ref. JTst by examining both the rotational symmetry of the static quark potential and the scaling behavior of light 
hadron mass ratios. 

Preliminary results for the topological susceptibility based on a first analysis at our intermediate lattice spacing 
have been published in Ref. In this article we present the final analysis and results at all gauge couplings. 

The identification of dynamical quark effects requires a comparison with pure SU(3) gauge theory where sea quarks 
are absent. We therefore supplement our study of topology in full QCD by a set of simulations of SU(3) gauge theory 
with the same RG-improved gluon action at a similar range of lattice spacings. 

The outline of this article is as follows. In Sec. |l| we give details on numerical simulations and measurements of the 



topological charge. Results for the topological susceptibility are presented in Sec. HI where we discuss the continuum 



extrapolation in pure gauge theory, as well as the quark mass dependence in full QCD. Conclusions are summarized 



in Sec. IV 



II. COMPUTATIONAL DETAILS 



A. Gauge configurations 

Gauge configurations incorporating two degenerate flavors of dynamical nuarks have been generated by the CP- 
PACS full QCD project. For gluons we employed an RG-improved action ||l^ of the form, 

^RG = ^ I 3.648 K^\^) - 0-331 E K^^W ' (2) 

where M^^^^ and W^^^ are the plaquette and rectangular Wilson loop. For the quark part we adopted the clover 
quark action with a mean fleld improved clover coefficient csvv = P^^^^, and the plaquette P calculated in 
perturbation theory at one loop P — 1 — 0.8412/3^^. This choice is based on the observation that measured values of 
the plaquette (P) are well approximated by the one-loop estimate and that csw determined in this way is close 
to its one- loop value [ p^ . 
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Three sets of gauge configurations have been generated at bare gauge couplings /? = 1.8, 1.95 and 2.1 corresponding 
to the lattice spacings a « 0.22, 0.16 and 0.11 fm. Lattices of size x T = 12^ x 24, 16^ x 32 and 24^ x 48 have been 
used, for which the physical lattice size remains approximately constant at La w 2.5 fm. At each (3, runs are carried 
out at four values of the hopping parameter k chosen such that the mass ratio of pseudoscalar to vector mesons takes 
TOps/mv « 0.8, 0.75, 0.7 and 0.6. 

In Table | we give an overview of the parameters and statistics of the full QCD runs. Technical details concerning 
the configuration generation with the Hybrid Monte Carlo (HMC) algorithm and results for the light hadron spectrum 
are presented in Ref. |lj]. Runs were made with a length of 4000-7000 HMC unit-trajectories per sea quark mass. 
Topology measurements are made on configurations separated by 10 HMC trajectories at /? = 1.8 and 1.95 and by 
5 trajectories at /? = 2.1. The number of measurements A^Moas and the separations iVskip are listed for each run in 
Table |. 

We supplement the study of topology in full QCD by simulations of pure SU(3) gauge theory with the RG-improvcd 
action of Eq. (|^). Configurations are generated at 5 values of (3 with lattice spacings 0.09 fm < a < 0.27 fm as listed in 
Table |l[ For the three larger gauge couplings lattices of size 8^, 12"* and 16'' are used so that the physical lattice size 
remains approximately constant at La ~ 1.5 fm. While this is smaller than the sizes in the full QCD runs, it has been 
a standard size employed in recent studies of topology in SU(3) gauge theory [p0|-p3|. It has also been shown |2|] 
that the instanton size distribution does not suffer from significant finite volume effects on a lattice of this size. For 
the two smaller gauge couplings we keep lattices of size 8^. Simulations are carried out with a combination of the 
pseudo-heat-bath algorithm and the over-relaxation algorithm mixed in a ratio 1 : 4. For each /3 we create 500-2000 
independent configurations separated by 100 iterations. 



B. Topological charge operator 



The topological charge density in the continuum is defined by 

1 



Q{x) 



327r2 



(3) 



and the total topological charge Q is an integer defined by the integrated form 

Q = Jd^x Q{x). (4) 
On the lattice we use the field theoretic transcription of this operator which has the standard form 

with the lattice charge density defined by 



1 



.,,Tr(C;,(x)C;,(x)) 



32^2^'^- 

In this expression the field strength on the lattice is defined through the clover leaf 



(6) 



(7) 



An improved charge operator can be constructed by additionally calculating a rectangular clover leaf made out of 
1x2 Wilson loops. 



and combine them to the charge density 



327r 



ypa Tr [C^^{x)Cf„{x)) 



(8) 



(9) 
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The improved global charge is then defined through 

Qimp = {coQLi^n) + ciQf (a;„)} . (10) 

The standard charge operator of Eq. (H) has O(a^) discretization errors. With the choice cq — 5/3 and 
ci = —1/12 |2^j2^ in Eq. ( p^ the leading order terms are removed for classical instanton configurations and 
discretization errors become 0(a^). 



C. Cooling 



The topological charge operators of Eq. (||) or (^ are dominated by local fluctuations of gauge fields when measured 



removes the 
ying long-range 



on thermalized lattice configurations and their value is generally noninteger. The cooling method 
ultraviolet fluctuations by minimizing the action locally while not significantly disturbing the underl 
topological structure. 

In full QCD one might consider cooling with the full action including the fermionic part. We refrain from this 
because it would lead to solutions of the classical equations of motion of the effective action, obtained by integrating 
out fermion fields pTj . These are different from instantons which are solutions of the classical equations of motion of 
the gauge action only. Moreover, cooling would become a non-local process. 

In principle any lattice discretization of the continuum gauge action can be used for smoothing gauge configurations 
by cooling. However, lattice actions generally do not have scale invariant instanton solutions. The standard Wilson 
plaquette action discretization of a continuum instanton solution with radius p, for example, behaves for a <ti p <ti L 
as 

^piaq = 5cont |l - ^(a/p)' + O ((a/p)4)| . (11) 

Under cooling with the plaquette action, instantons therefore shrink and disappear when the cooling is applied too 
long. To improve on this we use for cooling a gluon action of the generic form 



Co 



^ w;,^\x)+c,j2wi:^'ix)Y (12) 



where the coefficients cq and ci satisfy the normalization condition cq + 8ci — 1. We employ the two choices 

Co = 5/3, ci = -1/12 LW action, (13) 

and 

Co — 3.648, ci = —0.331 RG-improved action. (14) 

The tree-level improved Symanzik action by Liischer and Weisz [QUI of Eq. ( |l3|) has reduced breaking of instanton 
scale invariance given by [ p8[ , 

^LW = ^cont |l - + O ((a/p)6)| , (15) 

while still not admitting stable instantons under cooling. For the RG-improved action of Eq. ( p^ the sign of the 
leading order term is changed |28| , 

Skg = ^cont |l + ^(a/p)' + O ((a/p)4)| . (16) 

The flip of the sign leads to a local minimum of the action where stable lattice instantons can exist |2^ . 

Cooling with the RG-improved action or the LW action can lead to different values of the topological charge since 
instantons with a radius of the order of the lattice spacing can be either destroyed or stablized. The ambiguity is only 
expected to vanish when the lattice is fine enough. We test this explicitly by using both actions for cooling and treat 
differences as a systematic error of the cooling method. 



4 



A cooling step consists of the minimization of the local action for three SU(2) subgroups at every link of the lattice 
using the pseudo-heat-bath algorithm with [3 = oo. We have made 50 cooling steps for every configuration, measuring 
the topological charge after each step. 

We have investigated the deviations from integer topological charge as a function of the number of cooling steps, 
the topological charge operator, and the coupling constant for our simulations of pure SU(3) gauge theory. In Fig. |l| 
we show the distribution of the topological charge at the intermediate gauge coupling oi j3 — 2.227. The distribution 
is peaked at quantized but noninteger values of Q. The peaks are already well separated after 10 cooling steps and 
the widths of peaks further decrease with increasing number of cooling steps. At the same number of cooling steps, 
peaks are narrower for the improved charge operator Qimp than for the naive form Qst ■ Centers of peaks are located 



below integer values. Cooling and improvement of charge operator move them closer to integers. In Table III we 
list the ratio between center of peaks and integer charge, found to be independent of Q, for all gauge couplings and 
after 10, 20, or 50 cooling steps with the RG-improved action. The ratio moves closer to unity with increasing gauge 
coupling, increasing number of cooling steps and when the charge operator is improved, showing that the difference 
from integer is a finite lattice spacing effect. After 20 cooling steps, centers of peaks of Qimp do not differ from integer 
by more than 6% even at the coarsest lattice spacing. Because of its superiority we only use Qimp, rounded to nearest 
integer, in the following. 

In Fig. 1^ we plot {Q"^) measured in pure SU(3) gauge theory as a function of the number of cooling steps for the 
two cooling actions. Cooling with the two actions leads to quite different values of (Q^) at coarser lattice spacings. 
The difference decreases with increasing coupling constant and almost vanishes on the finest lattice. 



We quantify the difference between cooling with the two actions by calculating the linear correlation coefficient 



( nRG _ nRG ^ / qim _ qi 

I ^!imp ^imp J I ^imp 



(17) 



after 10, 20, or 50 cooling steps. For the evaluation of Eq. (|17|) we substitute charges before rounding to integers 




'3imp ) ) ( \Qimp Qimp 



Values of r are listed in Table IV. The correlation between topological charge after cooling with the RG-improved 
or the LW action decreases with increasing number of cooling steps. Even at the coarsest lattice spacing and after 
50 cooling steps, however, there is a strong correlation with r = 0.84. With decreasing lattice spacing r approaches 
unity and charges are highly correlated on the finest lattice. These features agree with our naive expectations. 

Since (Q^) has an approximate plateau after 20 cooling steps we use this as central value. (Q^) is listed for pure 
SU(3) gauge theory in Table and for full QCD in Table 0. The first quoted error is statistical. The second error 
expresses the uncertainty of choosing the number of cooling steps by taking the largest difference between (Q^) after 
20 cooling steps and after more cooling steps up to 50. 

D. Full QCD time histories 

Decorrelation of topology is an important issue in the simulation of full QCD since the topological charge is one 
of the quantities which is expected to have the longest autocorrelation with the HMC algorithm. In simulations with 
the Kogut-Susskind quark action it was found that topological modes have a very long autocorrelation time . 

In Figs. |[ ^ and || we plot time histories of Qimp after 20 cooling steps calculated for our full QCD runs at all sea 
quark masses. Autocorrelation times are visibly small even at the smallest quark masses. For (3 = 1.80 and (3 = 1.95 
the topological charges measured on configurations separated by 10 HMC trajectories are well decorrelated, and hence 
the integrated autocorrelation time is smaller than 10 trajectories. Correspondingly, errors are independent of the bin 
size when employing the binning method. At /3 = 2.10, where the charge is measured at every fifth HMC trajectory, 
we find integrated autocorrelation times of 5-6 configurations, corresponding to 25-30 HMC trajectories. This is 
comparable to, but somewhat smaller than, recent results reported for the Wilson |^ or the clover quark action 
For error estimates throughout this paper we use bins of 10 configurations, corresponding to 50 HMC trajectories, at 
/3 = 2.10 and no binning for the two other couplings. 

A related issue is the ergodicity of HMC simulations. In Figs. |^, ^ and ^ we show histograms of the topological 
charge. They are reasonably symmetric around zero and the distribution can be approximately described by a 
Gaussian, also plotted in the figures. Ensemble averages (Q), listed in Table 0, are consistent with zero or deviate 
at most three standard deviations of statistical error at /3 = 2.1 and k = 0.1374. We conclude that topology is well 
sampled in our runs. 
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E. Scale determination 



To fix the scale we use the string tension a or the Sonimer parameter tq of the static quark potential. Full 
QCD values of ro have been determined in Ref. |jl^ and are reproduced in Table |. The analysis of the static quark 
potential in pure SU(3) gauge theory of this work parallels the one in Ref. [^. We list a and ro in Table ||. The 
dependence of the dimensionless string tension \foa on the gauge coupling is shown for pure SU(3) gauge theory in 
Fig. ^ together with previous results ofRefs. ||l|,||,|4). Data are consistent with previous determinations, and extend 
the domain of results to smaller values of /?. 

We fit the string tension data of Fig. ^ using an ansatz proposed by AUton jsst , 

- {1 + C2a{p,f + C4a(/3)^} /cq, a{(5) = /(/?)//(/? = 2.4), (18) 
where /(/3) is the two-loop scaling function of SU(3) gauge theory, 

We obtain the best fit at 

Co = 0.5443(97), C2 = 0.390(38), C4 = 0.049(12), (20) 
with good /^DF = 19.3/19. The fit curve plotted in Fig. ^ reproduces the data very well. 



III. TOPOLOGICAL SUSCEPTIBILITY 
A. Pure SU(3) gauge theory 

The topological susceptibility 

..^^ m 

in pure SU(3) gauge theory is converted to the dimensionless number xt^Q using measured values of the Sommer scale 
ro and is quoted in Table Statistical errors of {Q^) and ro and the systematic error related to the choice of the 
number of cooling steps are added in quadrature. 

We plot XtTo as a function of a^/r§ in Fig. 0. Results obtained with the two cooling actions are significantly different 
from each other at coarser lattice spacings. As expected, they move closer together toward the continuum limit. On 
the finest lattice the difference almost vanishes. Since data exhibit a curvature, we attempt continuum extrapolations 
including the leading scaling violation term of 0{a^) and the next higher order term of 0{a^). We obtain 

4 _ J 0.0570(43) + 0.049(61)aVr2 - 0.44(19)a'^/r4 cool with Srg ,r,r,. 
^ \ 0.0602(43) - 0.072(68)aVr2 - 0.19{22)a'^ /rf^ cool with Sun, ^ ' 

with x^ /Njjp — 2.2 and 1.4 respectively. Fit curves plotted in Fig. ^ follow the data well. In the continuum limit 
Xtro obtained with the two cooling actions differ by about one standard deviation of statistics. 

In Fig. we also plot Xt normalized by the string tension. Data behave similar to the one normalized by ro. A 
continuum extrapolation of the same form as above leads to 

, 2 _ / 0.0333(27) + 0.004(29)cra2 - 0.103(67)cr2a4 cool with S'rg ,„o^ 
^ \ 0.0347(27) - 0.040(31)CTa2 - 0.041(76)(72a4 cool with S^lw, ^ ' 

with x^ lN£)p — 1.5 and 0.8 respectively. 

To set the scale we use ro — 0.49(3) fm or %/a — 440(30) MeV where the errors in parentheses are our estimates of 
uncertainty of these quantities which are not directly measurable in experiments. Employing Xtrg from cooling with 
the RG-improved action as the central value, we obtain for the topological susceptibility in pure SU(3) gauge theory, 

xJ/* = 197(4)(t;5)(tO)(12) MeV, (24) 

where the first error is statistical, the second is associated with the uncertainty from the cooling action, the third 
reflects the difference from using ro or to set the scale, and the last comes from the uncertainty in ro. 

Our value of x^^ is in good agreement with recent determinations by several groups using different methods pO|-|2^] 
as well as with the Witten-Veneziano relation [||], xt = jl{r(i%, + m\ - 2m\)l2Nj k, (180 MeV)". 
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B. Full QCD 



Topological susceptibilities obtained in fullQCD runs and normalized by tq measured for the same sea quark mass 
are collected in Table VI. In Figs. ||, || and ^ they are plotted as a function of (mps'^o)'^- As in pure SU(3) gauge 



theory, data obtained with the two cooling actions differ from each other at /3 = 1.8 where the lattice is coarsest but 
are consistent with each other within error bars at /3 = 2.1. The quark mass dependence is similar between the two 
cooling actions at all the /? values. 

For comparison we also plot in Figs. ^, ^ and |l^ susceptibilities in pure SU(3) gauge theory obtained by cooling 
with the RG-improved action. In full QCD vq changes together with the sea quark mass or mps in Table |[ The 
topological susceptibility in pure gauge theory is a decreasing function of a^/r§, and the value corresponding to full 
QCD at the same vq is therefore not constant when mpg changes. We take this into account by using the interpolation 
formula of Eq. ( p^ ) and the linear fit of l/ro as a function of rripg in Ref. and calculate Xt^a at matching values 
of 7'o. We arrive at the one standard deviation error band of the susceptibility in pure gauge theory plotted as the 
light shaded area in Figs. H, ^ and ^ An increasing tendency with decreasing quark mass is manifest at /3 = 1.8, 
whereas at (3 = 2.1 the shaded error band is very flat. 

The topological susceptibility in full QCD is consistent with that of pure gauge theory at the heaviest quark mass 
for P = 1.8 and 1.95, but smaller by two standard deviations for f3 = 2.1. Values at intermediate quark masses are 
consistent or slightly smaller. At the smallest quark mass the topological susceptibility in full QCD is suppressed 
compared to the pure gauge value. The decrease is, however, contained within 15% or one to two standard deviations 
at /? = 1.8 and 1.95, which is marginal. A clearer decrease by 41%, corresponding to seven standard deviations, is 
observed at /3 = 2.1. 

We investigate if the small suppression due to dynamical quarks at the two coarser lattice spacings is against 
expectations by comparing the behavior of xt with the prediction of Eq. (|l|) for vanishing quark mass. Using the 
Gell-Mann-Oakes-Renner relation ||37| , |38|| 

S = ^, (25) 
with normalized to be 132 MeV in experiment, Eq. (Q) can be rewritten as, 

X..o^ = ^^^^^^^^^^^ (26) 

In Ref. ||l^ pseudoscalar decay constants /ps and Sommer scale tq have been determined for all gauge couplings and 
fitted as functions of mpg. Using the fits /pg(rnpg) and »'o(™ps) we calculate the one standard deviation error band 
of Eq. (|2^) and plot it as dark shaded area starting at zero in Figs. |[ ^ and |l^. We plot the same prediction evaluated 
with measured values of /ps and at physical quark masses as dotted line. Differences between the band and the line 
are of order mpg. Sizable scaling violations in /pg have been observed in Ref. with /^r — 195(5) MeV (/3 — 1.8), 
157(7) MeV (/3 — 1.95) and 131(7) MeV (/3 — 2.1) if the scale is determined by the p meson mass. Correspondingly, 
the slope of the prediction Eq. ( p6| ) shows a variation with (3. 

The susceptibility in full QCD at the smallest quark mass lie between the shaded band and the dotted line of 
Eq. (p^). Interestingly, the smallest simulated quark masses at /3 = 1.8 and 1.95 lie roughly in the region where the 
small mass prediction and pure SU(3) gauge theory cross. A stronger suppression of the topological susceptibility at 
f3 = 2.1, on the other hand, occurs at a quark mass somewhat below the crossing point. This may be an indication 
that the runs at /3 = 2.1 reach quark masses where a suppression compared to pure SU(3) gauge theory can be 
expected. The exact location of the cross-over region depends, however, on the magnitude of higher order terms in 
Eq. ( p6| ) and the lattice value of fps- Simulations at lighter quark masses will therefore be helpful to clarify whether 
the interpretation described here is correct. 



IV. DISCUSSION AND CONCLUSIONS 



We have studied the topological susceptibility as a function of quark mass and lattice spacing in two-flavor full 
QCD using a field theoretic definition of the topological charge together with cooling. 

We have shown that an improved charge discretization can be defined which produces charges close to integers. 
The stability of lattice instantons differs between two actions used for cooling, which leads to different values of the 
topological charge at coarse lattice spacings. We have confirmed that the difference decreases with decreasing lattice 
spacing and vanishes in the continuum limit. Our investigation of time histories of the topological charge in full QCD 
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have shown that autocorrelations are reasonably short and that our runs are long enough to sample topology well. 
These analyses support our belief that systematic errors of the cooling method are kept under control, and that our 
lattice measurements indeed reflect topological properties of the QCD vacuum. 

The quark mass dependence of the topological susceptibility Xti^o iii full QCD is found to be flat or even increase 
with decreasing quark mass at (3 — 1.8 and 1.95, and a clear decrease is only observed aX (3 — 2.1. A comparison with 
pure gauge theory at corresponding rg shows that Xt^o f^l^ QCD is consistent with pure gauge theory at heavier 
quark masses but suppressed at the lightest quark mass of our simulation. At the same time, the susceptibility at 
the lightest quark masses are in agreement with the prediction of the anomalous Ward-Takahashi identity for U(l) 
chiral symmetry for small quark masses when lattice values for the pseudoscalar decay constant are employed. These 
results suggest that our lightest simulated quark masses lie around the transition region where a suppression due to 
sea quarks is expected to set in. 

Recently several alternative theoretical explanations have been suggested as to why the topological susceptibility 
in lattice full QCD might appear less suppressed than expected for small quark masses. It has been pointed out |3^] 
that a large enough volume with VT^niq » 1 Q is necessary for Eq. (0) to be valid. Since we employ a large lattice 
size of La « 2.5 fm and quark masses with >40 MeV this condition is always fulfilled. It has also been argued 
that subtleties exist in the flavor singlet lattice Ward-Takahashi identities when the Wilson or clover fermion action 
is employed so that counter-terms are needed for the correct chiral behavior of the topological susceptibility [Q. 
Simulations at lighter sea quark masses and further theoretical analyses are needed to examine whether such an 
explanation is required for understanding the quark mass dependence of Xt in full QCD. 
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TABLE I. Overview of full QCD simulations. The lattice spacing a is fixed by the vector meson mass at the physical quark 
mass and Mp = 768.4 MeV. 





L-^ X T 


csw 


a [fm] 


La [fm] 


K 


mpso 


mps/mv 


ro/a 


N-Meas 




A'^Bin 


1.80 


12^x24 


1.60 


0.2150(22) 


2.580(26) 


0.1409 


1.15601(61) 


0.807(1) 


1.716(35) 


650 


10 


1 












0.1430 


0.98267(89) 


0.753(1) 


1.799(13) 


522 


10 


1 












0.1445 


0.82249(82) 


0.694(2) 


1.897(30) 


729 


10 


1 












0.1464 


0.5306(17) 


0.547(4) 


2.064(38) 


409 


10 


1 


1.95 


16^x32 


1.53 


0.1555(17) 


2.489(27) 


0.1375 


0.89400(52) 


0.804(1) 


2.497(54) 


681 


10 


1 












0.1390 


0.72857(68) 


0.752(1) 


2.651(42) 


690 


10 


1 












0.1400 


0.59580(69) 


0.690(1) 


2.821(29) 


689 


10 


1 












0.1410 


0.42700(98) 


0.582(3) 


3.014(33) 


488 


10 


1 


2.10 


24^x48 


1.47 


0.1076(13) 


2.583(31) 


0.1357 


0.63010(61) 


0.806(1) 


3.843(16) 


800 


5 


10 












0.1367 


0.51671(67) 


0.755(2) 


4.072(15) 


788 


5 


10 












0.1374 


0.42401(46) 


0.691(3) 


4.236(14) 


779 


5 


10 












0.1382 


0.29459(85) 


0.576(3) 


4.485(12) 


789 


5 
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TABLE IL Overview of pure SU(3) simulations. The lattice spacing a is determined using ^/a = 440 MeV. Numbers in 
brackets for A^conf indicate the number of configurations used for potential measurement. 



/3 


L '-^ X T 


a [fm] 


La [fm] 




ro/a 




2.047 


8^x8 


0.2726(19) 


2.181(15) 


0.3695(52) 


1.8978(59) 


500 


2.110 


8^x8 


0.2439(10) 


1.951(8) 


0.2958(24) 


2.1399(53) 


1000 


2.227 


8^x8 


0.1905(10) 


1.524(8) 


0.1805(19) 


2.738(11) 


2000 


2.461 


12^x12 


0.1259(7) 


1.511(9) 


0.07885(90) 


4.089(14) 


900 


2.659 


16^x16 


0.0931(9) 


1.489(14) 


0.04311(84) 


5.556(30) 


700(495) 



TABLE III. Ratio between center of peak of the topological charge distribution and integer charge after 10/20/50 cooling 

steps with the RG-improved action. At missing numbers no clearly separated peak structure could be identified. 

(3 standard Q improved Q 

2.047 — /0.77/0.85 — /0.94/0.97 

2.110 — /0.80/0.87 0.89/0.94/0.98 

2.227 0.78/0.83/0.88 0.94/0.96/0.98 

2.461 0.85/0.89/0.92 0.97/0.98/0.99 

2.659 0.89/0.92/0.94 0.98/0.99/0.995 



TABLE IV. Correlation coefficient r between Qimp obtained after 10, 20 or 50 cooling steps with the RG-improved or the 

LW action. 

(3 10 steps 20 steps 50 steps 

2.047 0.90(1) 0.86(1) 0.84(1) 

2.110 0.923(5) 0.886(7) 0.871(8) 

2.227 0.961(2) 0.942(3) 0.931(3) 

2.461 0.991(1) 0.986(2) 0.982(3) 

2.659 0.9982(5) 0.9978(7) 0.9970(9) 
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TABLE V. Topological susceptibility in pure SU(3) gauge theory. For (Q^) the first error is statistical, and the second error 
is an estimate of systematic error related to the choice of the number of cooling steps. For Xf'o the two errors and the statistical 
error of ro are added in quadrature. 



/3 


cool with RG-improved action: 

(Q) (Q'} 


Xtro 


(Q) 


cool with LW action: 


Xtrf) 


2.047 
2.110 
2.227 
2.461 

2.659 


0.05(15) 
0.121(93) 
-0.043(46) 
0.139(68) 
0.067(76) 


12.07(69)(+72) 
8.61(39)(+13) 
4.24(13)(0) 
4.12(21)(0) 
4.08(22)(0) 


0.03821^2 
0.044lt^; 

0.05821^0 

0.0555^29 
0.0593^34 


0.00(13) 
0.054(82) 
-0.042(43) 
0.123(66) 
0.073(76) 


8.50(52)(-80) 
6.74(31)(-38) 
3.71(12)(-18) 
3.93(20) (-4) 
4.06(22)(-l) 


0.0269^30 
0.0345125 
0.050913? 
O.O530I28 
O.O590I34 



TABLE VL Topological susceptibility in full QCD. Meaning of errors are the same as in Table 









cool with RG-improved action: 




cool with LW action: 




/3 




K 


(0) 




Xtrf) 


(Q) 




Xtro 


1.80 





1409 


0.41(44) 


123.4(6.5)(+11.9) 


0.02581^^ 


-0.04(37) 


87.9(4.6)(-3.3) 


0.0184li^ 







1430 


0.10(49) 


125.1(7.8)(+10.5) 


O.O3I6I22 


0.20(40) 


85.2(5.5)(-3.2) 


O.O215II7 







1445 


0.46(40) 


119.3(6.0)(+10.7) 


0.037113^ 


0.53(33) 


77.7(4.2)(-3.0) 


O.O241I22 







1464 


0.18(46) 


85.0(5.9)(+6.8) 


0.0372l3^ 


0.20(38) 


58.1(3.9)(-0.7) 


O.O254I25 


1.95 





1375 


1.09(52) 


186.4(9.7)(+10.7) 


0.0553l^g 


1.35(46) 


148.5(8.4)(+1.6) 


0.0440l4^ 







1390 


0.42(43) 


127.0(6.5)(+9.3) 


O.O478I39 


0.27(39) 


104.2(5.7)(-0.6) 


O.O393I33 







1400 


-0.33(39) 


106.3(5.7)(+8.3) 


0.05141^^ 


-0.02(34) 


78.4(4.4)(+1.9) 


O.O379I2? 







1410 


0.61(40) 


76.5(4.7)(+5.0) 


O.O482I3® 


0.27(37) 


65.2(3.9)(+1.4) 


0.04111^1 


2.10 





1357 


0.96(88) 


146.4(11.0)(+6.8) 


O.O48II37 


1.17(88) 


137.4(10.6)(+2.6) 


0.0452l^g 







1367 


-0.5(1.0) 


150.7(16.6)(+5.7) 


0.06241^^ 


-0.6(1.0) 


137.6(15.6)(+3.4) 


0.0570165 







1374 


-2.52(81) 


102.2(9.7)(+3.3) 


0.0496l4° 


-2.61(82) 


98.0(9.7)(+0.9) 


0.04721^2 







1382 


-0.29(63) 


56.5(5.6)(+0.9) 


0.0344li^^ 


-0.33(61) 


52.6(5.0)(+0.5) 


0.03211^1 
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FIG. 1. Topological charge distribution at /3 = 2.227 after various numbers of cooling steps with the RG-improved action 
and for two definitions of the topological charge. Figure a) shows the whole distribution while Fig. b) is an enlargement of the 
first three peaks. 
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FIG. 2. Expectation value of the topological charge squared as a function of the number of cooling steps for two different 
cooling actions. 
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FIG. 4. Time histories and histograms in fuU QCD at /3 = 1.95. 
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FIG. 5. Time histories and histograms in full QCD at /3 = 2.10. 
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FIG. 6. String tension in pure SU(3) gauge theory as a function of the gauge coupling. Circles represent data from 
Refs. p^|3^ , p^ while squares are obtained in the present work. The solid line represents a fit with Eq. (p^. 
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FIG. 7. Continuum extrapolation of the topological susceptibility in pure SU(3) gauge theory. 
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FIG. 8. Topological susceptibility in full QCD at /3 = 1.80. The light shaded region indicates the one standard deviation 
error band for pure SU(3) gauge theory, cooled with the RG-improved action, at corresponding values of ro. The darker shaded 
region starting at zero is the one standard deviation error band of the small mass prediction of Eq. (^) evaluated with measured 
values of /ps(?nps) and r-o(mps) while the dotted line is the same prediction evaluated with measured values of /ps and ro at 
physical quark masses. 
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FIG. 9. Topological susceptibility in full QCD at /3 = 1.95. Symbols are the same as in Fig 
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FIG. 10. Topological susceptibility in full QCD at /9 = 2.10. Symbols are the same as in Fig. 
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